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This  brief  DEPSCoR/AFOSR  sponsored  research  grant  has  been  used  to 
continue  the  extension  of  our  research  effort  on  the  asymptotic  theory  of 
ultrawideband  electromagnetic  pulse  propagation  through  temporally 
dispersive,  attenuative  media  to  spatially  inhomogeneous,  conducting, 
gyrotropic  media,  with  particular  emphasis  on  pulsed  wave  propagation 
through  the  ionosphere  and  magnetosphere.  This  funded  research  has 
provided  a  detailed  asymptotic  description  of  ultrawideband 
electromagnetic  pulse  propagation  through  a  semiconducting  medium  as 
described  by  the  Debye  model  with  static  conductivity.  This  analysis 
shows  that  the  presence  of  conductivity  does  not  eliminate  the 
appearance  of  the  Brillouin  precursor  that  is  associated  with 
ultrawideband  pulse  propagation  through  a  Debye  model  dielectric  and 
dominates  the  propagated  field  evolution  as  the  propagation  distance  z 
exceeds  a  single  absorption  depth  ^  at  the  input  pulse  carrier  frequency, 

decaying  only  as  z*  as  z^>  With  conductivity  present,  the  Brillouin 
precursor  now  experiences  some  exponential  attenuation  in  addition  to 
its  characteristic  z *  algebraic  decay.  Nevertheless,  our  analysis  has 
shown  that  this  exponential  attenuation  is  less  than  that  experienced  by 
the  remainder  of  the  pulse.  As  a  consequence,  the  Brillouin  precursor 
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remains  as  the  pulse  component  that  experiences  minimal  attenuation  in 
a  dispersive  attenuative  medium,  even  in  the  presence  of  static 
conductivity. 

The  results  of  this  funded  research  have  been  published  in  the 
Proceedings  of  the  Fourth  LASTED  International  Conference  on 
Antennas,  Radar,  and  Propagation  that  was  held  in  Montreal.  Canada  in 
May  2007.  A  copy  of  this  paper  is  attached  to  this  report. 
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ULTRAWI DEBAND  PULSE  PENETRATION  IN  A  DEBYE  MEDIUM  WITH 

STATIC  CONDUCTIVITY 
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ABSTRACT 

The  high  frequency,  low  power  characteristics  of  ultra- 
wideband/short-pulse  signals  exhibit  uniquely  promising 
features  with  applications  to  remote  sensing  of  terrestrial 
objects  from  satellites,  foliage-penetrating  radar,  as  well 
as  the  study  of  biological  tissue  exposed  to  ultrawideband 
pulses.  Indeed,  previous  studies  of  ultrawideband  elec¬ 
tromagnetic  pulse  propagation  through  dispersive,  noncon¬ 
ducting  media  has  shown  the  existence  of  a  so-called  Bril- 
iouin  precursor  whose  peak  amplitude  only  decays  alge¬ 
braically  with  propagation  distance.  However,  materials 
such  as  the  ionosphere,  foliage  and  biological  tissue  ex¬ 
hibit  conductivity.  In  this  paper,  we  show  that  a  Dcbye- 
model  material  with  static  conductivity  does  indeed  support 
a  Brillouin  precursor,  but  that  this  precursor  now  attenu¬ 
ates  exponentially  with  propagation  distance  and  not  just 
algebraically.  Nevertheless,  we  show  that  it  is  still  advan¬ 
tageous  to  track  the  Brillouin  precursor  in  remote  sensing 
applications  because  its  attenuation  is  less  than  the  attenu¬ 
ation  of  the  main  signal. 

KEY  WORDS 

wave  penetration,  ultrawideband  pulse  propagation,  pre¬ 
cursors 

1  Introduction 

An  ultrawideband  (UWB)  electromagnetic  pulse  is  a  pulse 
whose  bandwidth  spans  tens  of  megahertz  to  several  gi¬ 
gahertz.  In  the  time  domain,  this  corresponds  to  signals 
that  have  a  pulse  duration  on  the  order  of  1 00  picoseconds 
(ps)  or  less  (short- pulse),  or  signals  that  have  a  sufficiently 
rapid  turn-on  lime  also  of  order  100  ps  or  less.  Ullrawide- 
band  technologies  are  being  touted  by  some  as  one  of  the 
most  promising  new  technologies  of  our  time  because  their 
high  frequency,  low  power  characteristics  offer  security,  in¬ 
creased  penetration  capabilities,  immunity  from  multipath 
interference,  precise  positioning  capability  and  high  data 
rates.  Advances  in  generating  such  pulse  types  will  en¬ 
able  ultrawideband  pulses  to  be  used  in  an  increasing  num¬ 
ber  of  leading-edge  technologies  such  as  remote  sensing  of 
terrestrial  objects  from  satellites,  foliage- penetrating  radar, 
next-generation  cellular  telephones,  as  well  as  the  study  of 
biological  tissue  exposed  to  ultrawideband  pulses.  All  ap¬ 
plications  of  UWB  will  involve  the  propagation  of  UWB 


pulses  through  dispersive  material  such  as  the  earth's  sur¬ 
face,  trees,  walls  or  the  ionosphere.  Previous  studies  ( l  ]  of 
UWB/ short-pulse  electromagnetic  signals  through  disper¬ 
sive,  nonconducting  material  has  show  n  the  existence  of  a 
Brillouin  precursor  wfhose  peak  amplitude  decays  only  al¬ 
gebraically  with  propagation  distance;  a  considerable  at¬ 
tribute  for  application  purposes.  However,  many  physi¬ 
cal  materials,  such  as  the  ionosphere,  foliage  and  biolog¬ 
ical  tissue,  arc  conductive.  The  dielectric  permittivity  of 
a  conductive  material  possesses  a  simple  pole  at  the  ori¬ 
gin  u;  —  0,  which  is  the  very  point  that  provides  the  peak 
amplitude  of  the  Brillouin  precursor  and  so  it  has  been  sug¬ 
gested  that  conductivity  prevents  the  formation  of  a  Bril- 
lioun  precursor.  In  this  paper,  we  present  an  analysis  of  an 
ultrawideband  electromagnetic  pulse  propagating  through 
a  Debye  material  [2]  with  static  conductivity  and  we  study 
the  effect  of  conductivity  on  the  formation  of  the  Brillouin 
precursor. 

2  Formulation 

Let  all  of  spaec  be  occupied  by  a  Debye-type  material  u  ith 
static  conductivity  whose  relative  complex  dielectric  per¬ 
mittivity  is  given  by 

At  .a 

*c(^)  —  ^  : - : - M-,  ( 1 ) 

I  —  VjJT  u; 

where  u)  represents  angular  frequency,  is  the  high  fre¬ 
quency  limit  of  the  dielectric  permittivity,  Ar  =  ts 
where  fs  =  r(f))  is  the  static  permittivity  of  Ihe  mate¬ 
rial,  r  is  the  characteristic  relaxation  time  of  the  mate¬ 
rial,  and  a  —  rro/fo  where  <7q  is  the  static  conductivity 
of  the  material  and  f0  is  the  permittivity  of  free  space. 
Here,  we  employ  ihe  values  =  15.24,  fs  —  23.10  and 
r  =  1,16  x  10-10s  which,  writh  *7q  zz  1  x  10~jomho/m  are 
appropriate  parameter  values  for  leafy  vegetation  [3]. 

Consider  a  linearly-polarized,  electromagnetic  plane- 
wave  traveling  in  the  positive  z -direct ion.  Let  the  electric 
field  on  the  plane  c  —  0  be  a  step- function  modulated  sinu¬ 
soid  of  fixed  carrier  frequency 

E(G,  i)  =  u(t)  sin(ijct}x,  (2) 

where  u(t)  is  the  Heaviside  step- function  and  u)c  is  the 
fixed  carrier  frequency  of  the  pulse.  For  the  examples  given 


here,  ivV  —  2tt  x  106  rad/s.  Then  the  electric  field  compo¬ 
nent  of  the  propagated  signal  on  any  plane  z  >  0  is  given 
exactly  by  the  Fourier- Lap  lace  integral  representation  [  I  j 


0) 


where  u  >  U  is  greater  that  the  abscissa  of  absolute  con¬ 
vergence  for  the  initial  electric  field  component  E{ 0?  f),  c 
is  the  speed  of  light  in  vacuum,  and  0(uj,  0)  is  the  complex 
phase  function  defined  as 


-  iw[n(u))  -  0],  (4) 

w  here  tj.(w)  -  is  the  complex  index  of  refrac¬ 

tion  and  0  =  ct/z  is  a  dimensionless  space-time  parameter. 
Here,  the  magnetic  permeability  is  always  taken  as  ft  =  L 
Note  that  in  the  integral  representation  (3),  D  —  uj*  +  iu,r 
is  complex. 

We  concern  ourselves  only  with  the  determination 
of  the  electric  field  component  because  the  magnetic  field 
component  of  the  transverse  plane-wave  can  be  determined 
by  the  relation  B  =  k(u/)  x  E,  where  k(u?)  =  (w/c)n(u;) z 
is  the  complex  wave  vector.  Note  that  our  initial  pulse  is  ul- 
irawideband  in  the  sense  that  its  spectrum  decays  as 
as  uj  —*  oo.  Unfortunately,  the  integral  (3)  has  no  known 
closed-form  solution  valid  for  all  values  of  0  (except  in  vac¬ 
uum).  Instead,  we  apply  uniform  asymptotic  techniques  to 
provide  an  analytic  approximation  to  the  propagated  field, 
the  accuracy  of  this  approximation  increasing  with  increas¬ 
ing  propagation  distance. 


J  Topology  of  the  Complex  Phase  Function 


On  this  branch,  0)}  is  negative,  as  characteristic  of 

a  lossy  medium. 

For  values  of  0  <  sftZ,  the  contour  of  integration 
appearing  in  the  integral  representation  of  the  propagated 
field  (3)  may  be  completed  in  the  upper  half  of  the  com¬ 
plex  tu-plane,  For  sufficiently  large  |u>],  the  real  part  of  the 
complex  phase  function  behaves  as 

*(w,0]  =  =  -*/*  (yfcZ-  0}  +  0(1),  (7) 

and  because  is  analytic  in  the  upper  half  plane  (ex¬ 

cluding  the  real  axis),  it  follows  from  residue  theory  that  the 
propagated  field  is  identically  zero  for  //  <  y/rZ-  One  may 
also  prove  that  the  propagated  field  is  zero  at  the  space-time 
point  0  =  yJTZ  by  rewritling  Eq.  (3)  as 

E(z,t)  =  [  j1  2 exp  [-0(hJ,0)1  do;  (8) 

and  then  using  a  similar  argumcnl  as  above.  Hence*  in  the 
following  discussion  wrc  consider  the  evaluation  of  Eq.  (3) 
for  values  of  0  >  y/TZ- 

The  saddle  points  of  the  complex  phase  function  are 
solutions  to  the  equation 

0)  =  i  [7i(lu)  -  0  4-  um'(u;)]  =  0,  (9) 

where  the  prime  denotes  differentiation  with  respect  to 
There  arc  four  solutions  to  (9),  but  one  solution  lies  along 
the  imaginary  axis  below  the  branch  point  and  may  thus 
be  ignored  since  it  is  not  needed  in  the  remaining  analysis. 
The  locations  of  the  other  three  saddle  points  may  be  ap¬ 
proximated  in  the  region  M  <  1  (r  for  <  0  <  oo  by 
the  expressions 


The  complex  phase  function  of  a  Debye  medium  with  static 
conductivity 
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(5) 


is  a  multi-valued  function  with  four  branch  points,  one  lo¬ 
cated  at  the  origin  and  the  remaining  located  at  u) 1,2.3  where 
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Note  that  the  radicand  appearing  in  (6)  is  always  positive 
for  physically  reasonable  parameter  values  £.«,  eoo.er  and  r. 
The  complex  phase  function  is  made  single-valued  by  con¬ 
struction  of  two  branch  cuts  along  the  imaginary  axis:  one 
connects  the  origin  and  wa,  the  other  connects  lU|  and  W3. 
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+  9  <  On-, 

v.,b)I/S.  »  =  <,0> 

\t{-‘-7TT’)’  e<fe; 

(A)''’.  «*««  <"> 

,-i*(7ife  +  1)'  *>*» 

"<0" 


(A)  •  *-“» 

I  +  f)  .  »  >  «o 


where 


D  —  ~  3too)  +  +  foo), 


(12) 


(13) 


and  0o  =  These  expressions  show  that  u7sp,.(0)  lies 
along  the  positive  imaginary  axis  for  all  0  >  y/ZZ  and 
approaches  the  branch  point  at  the  origin  as  0  -+  oo,  while 
u/5F,  a(0)  evolve  in  the  complex  lower  half  plane  for  < 
0  <  0q,  coalesce  into  a  second-order  saddle  point  along  the 


Figure  !  Evolution  of  the  saddle  points  of  0)  for  a 
Debye  medium  with  static  conductivity  <7^  =  l  x  !0_lD 
mho/m.  The  saddle  points  ^spi  l  3  (9)  are  depicted  by  the 
plus  sign,  cross  and  circle,  respectively. 


Figure  2.  The  Drillouin  precursor  of  a  step- function  mod¬ 
ulated  sine  wave  of  carrier  frequency  ujc  —  2tt  x  1()!‘  rad/s 
at  a  distance  of  into  a  Debye  medium  with  static  con¬ 
ductivity  on  =  1  x  10“ 10  mho/m  ♦ 


negative  imaginary  axis  between  the  branch  points  uja  and 
at  some  space-time  point  0S  Oq  ,  and  then  split  again 
into  two  first-order  saddle  points  with  wspt(0)  traveling 
down  the  imaginary  axis  toward  the  branch  point  as  9  — ^ 
oo  while  o;sf2(0)  travels  upward  toward  the  branch  point 
ui  as  0  -*  oo.  The  evolution  of  the  saddle  points  of  the 
complex  phase  function  for  a  Debye  medium  with  static 
conductivity  <r0  =  1  *  10” 10  mho/m  is  shown  in  Figure  1 . 
The  saddle  points  wspJi2,3  (0)  are  depicted  by  the  plus  sign, 
cross  and  circle,  respectively. 


4  Asymptotic  Approximation 

As  required  by  the  saddle-point  method,  the  contour  ap¬ 
pearing  in  Eq,  (3)  is  deformed  into  a  new  path  V(9)  such 
that  V{0)  is  continuous  in  6  and  £>(0)  passes  through  the 
dominant  (and  accessible)  saddle  points  of  *j>{u.\0)  while 
remaining  in  their  valleys.  For  the  analysis  here,  it  is  suffi¬ 
cient  to  consider  a  path  V{9)  that  passes  through  u j$px  (#) 
when  9  <  0O  and  ^sp3{0)  when  0  >  0O.  Of  course,  in 
the  deformation  of  the  original  Bromwich  contour  into  the 
path  P(0)t  the  branch  cut  extending  form  u>  -  0  to  is 
encircled  for  0  >  0O^  For  the  purpose  of  this  paper,  the 
integration  around  this  branch  cut  is  negligible  [4] 


4,1  Bril  ton  in  Precursor 

The  asymptotic  contribution  of  the  saddle  points  ljsfAQ) 
for  0  <  0O  and  ^sp2(0)  for  0  >  0O  yields  the  Brillouin 
precursor,  EB(z,t).  Use  of  the  saddle-point  method  to  the 


integral  representation  Eq.  (3)  yields 


27T!I?{wsP  -  LJC\] 


2  C7T 


■exp  ,  (14) 


where 


\«Sft(9),  for 8>0o. 


(15) 


An  example  of  the  Brilluoin  precursor  in  this  Debye 
medium  with  static  conductivity  00  =  1  x  10“  H>  mho/m  at 
a  propagation  distance  of  three  absorption  depths  is  given 
in  Figure  2,  Here,  z<i  is  one  absorption  depth  at  the  in¬ 
put  signal  carrier  frequency  uv  and  is  given  by  z<i  = 


4,2  Signal  Contribution 

The  signal  contribution  Ec(z ,  /)  to  the  propagated  field  so¬ 
lution  is  due  to  the  simple  pole  located  at  u>r  appearing 
in  (3),  When  S{$(u/c,0)}  ~  0)},  the  ef¬ 

fect  of  the  pole  on  the  saddle  point  u;sp(0)  must  be  taken 
into  account.  Also,  when  the  path  of  integration  P(0) 
crosses  the  pole,  the  appropriate  residue  contribution  must 
be  added  to  the  field,  A  direct  application  of  the  uniform 
theory  of  Felsen  [5]  and  Felscn  and  Marcuvitz  [6]  yields 
the  signal  contribution 

jvc(£,  £)  — 

-  |±i;rexp  0)j  erfc  ^T*A(0) 

+  ik/r 1 exp  }  .3  {*(*)}  §  0; 

(16) 


Figure  3,  The  signal  contribution  of  a  step- function  modu¬ 
lated  sine  wave  of  carrier  frequency  u-v  —  2 n  x  ]()6  rad/s 
at  a  distance  of  3ztt  into  a  Debye  medium  with  static  con¬ 
ductivity  cr0  =  1  x  10-10,  mho/m  , 


Figure  4.  Asymptotic  (dashed)  and  numerical  (solid)  solu¬ 
tions  to  the  propagated  electric  field  component  of  a  step- 
function  modulated  sine  wave  of  carrier  frequency  = 
2tt  x  It)6  rad/s  at  a  distance  of  3e(j  into  a  Debye  medium 
with  static  conductivity  =  1  x  10^ 10  mho/m. 


Er{z-  4)  — 

-  Tjf  |±i7rexp  ^0(wo0)]  erfc 

+W)^T'™P 

-  S  {zexp  )  ,0  {A(fl)}  =  l); 


(17) 


where 


A(fl)  =  ^(wspW.fi)  -0(^.9).  ( 1 8) 

with  erfc  {()  denoting  the  complementary  error  function, 
and  u)sp(b)  is  given  in  (15).  An  example  of  the  signal 
contribution  in  this  Debye  medium  with  static  conductiv¬ 
ity  (To  =  1  x  10” 10  mho/m  at  a  propagation  distance  of 
three  absorption  depths  is  given  in  Figure  3. 


4.3  The  Total  Field 

The  asymptotic  approximation  of  the  propagated  electric 
field  component  on  any  plane  z  >  0  is  given  by  the  sum 

E(zft)  =  EB(Ztt)  +  Ec(z:t).  (19) 

As  an  example  to  illustrate  the  accuracy  of  the  asymptotic 
approximation,  Figure  4  shows  both  the  asymptotic  solu¬ 
tion  and  a  numerical  solution  to  the  propagated  electric 
field  component  (3)  of  a  step- function  modulated  sine  wave 
of  applied  frequency  =  2tt  x  106  rad/s  traveling  through 
a  Debye  medium  with  static  conductivity  <r0  —  1  x  LQ“10 
mho/m  at  a  propagation  distance  of  3 zj.  The  two  results 
arc  almost  indistinguishable. 


5  Peak  Amplitudes 


The  penetration  capabilites  of  the  Brillouin  prccusor{l4) 
are  now  c%'aluatcd.  We  consider  three  different  values  for 
the  static  conductivity:  op  =  I  x  10“ ll,  1  x  1()^10T  1  x  lt)~^ 
mho/m*  For  each  value  of  conductivity,  the  peak  ampli¬ 
tude  of  the  Brillouin  precursor  is  recorded  at  the  relative 
distances  zjza  =  1  to  10  in  increments  of  one  and  from 
z/zd  =  10  to  100  in  increments  of  ten,  into  the  Debye 
material*  A  cubic  spline  is  fit  to  the  three  sets  of  data 
points  and  are  shown  in  Figure  5  by  the  dashed,  dotted 
and  dashed-dotted  curves,  respectively.  The  lower  solid 
curve  denotes  the  exponential  attenuation  exp)— z/'£d]  ex¬ 
perienced  by  the  carrier  frequency  of  the  pulse.  The  power 
of  the  observed  peak  amplitude  decay  is  found  by  taking 
the  logarithm  of  the  peak  amplitude  values  and  the  loga¬ 
rithm  of  the  relative  distances.  The  average  slope  of  these 
curves  is  calculated  and  then  plotted  in  Figure  6.  The 
average  slopes  for  the  three  levels  of  static  conductivity 
=  1  x  10"l\l  x  It)-10, 1  x  10_J1  mho/m  are  repre¬ 
sented  by  the  dashed,  dotted  and  dashed-dotted  curves,  re¬ 
spectively.  Evident  in  these  figures  is  the  rapid  decay  of  the 
carrier  frequency  of  the  pulse  and  the  minima!  impact  that 
the  static  conductivity  (at  least  for  these  values)  has  on  the 
peak  amplitude  values  of  the  Brillouin  precursor.  Hence, 
it  would  be  advantageous  to  track  the  Brillouin  precursor 
rather  than  the  carrier  frequency  in  remote  sensing  applica¬ 
tions  through  this  Debye  medium  w  ith  these  levels  of  static 
conductivity. 


6  Conclusion 


Figure  5.  Peak  amplitude  values  of  the  Brillouin  precusor 
plotted  as  a  function  of  distance  for  three  different  values 
for  the  static  conductivity:  cr0  —  1  x  10-U,  1  x  LG”10,1  x 
It}-9  mho/m  given  by  the  dashed,  dotted  and  dashed-dotted 
curves,  respectively.  The  lower  black  curve  denotes  the  ex¬ 
ponential  attenuation  exp[—  z/zd]  experienced  by  the  car¬ 
rier  frequency  of  the  pulse. 


Figure  6,  The  average  slopes  of  the  logarithm  of  the  data 
presented  in  Figure  5, 


An  asymptotic  approximation  to  the  propagated  electric 
field  component  of  a  step-function  modulated  sine  wave 
through  a  Debye  medium  with  static  conductivity  has  been 
presented  here.  The  accuracy  of  this  approximation  in¬ 
creases  with  increasing  propagation  distance.  The  asymp¬ 
totic  solution  consists  of  two  contributions:  the  Brillouin 
precursor  and  the  main  signal  contribution,  The  peak  am¬ 
plitude  points  of  these  contributions  are  plotted  as  functions 
of  relative  propagation  distance.  Both  components  have 
a  peak  amplitude  point  that  decays  exponentially,  but  the 
peak  amplitude  point  of  the  Brillouin  precursor  decays  less 
rapidly.  For  this  reason,  it  may  be  advantageous  to  track  the 
Brillouin  precursor  rather  than  the  carrier  frequency  in  re¬ 
mote  sensing  applications  through  this  Debye  medium  with 
these  levels  of  static  conductivity. 

7  Acknowledgements 

The  research  presented  in  this  paper  has  been  supported 
under  AFOSR  grant  #  000020247, 

References 

1 1]  K.  E.  Ough stun  and  G,  C  Sherman,  Pulse  Propagation 
in  Causal  Dielectrics,  Berlin:  Springer- Verlag,  1994. 

[2]  P,  Debye,  Polar  Molecules.  New  York:  The  Chemical 
Catalog  Company,  Inc.,  1929, 

[3]  R.  L.  Medina,  1  W.  Penn,  and  R.  A.  Albanesc,  "Di¬ 
electric  response  data  on  materials  of  military  conse¬ 
quence"  tech,  rep,,  Armstrong  Research  Laboratory, 
2001. 

[4]  For  the  material  parameters  used  here,  there  is  a  small 
theta  domain  in  which  the  path  of  integration  must 
pass  through  all  three  saddle  points  before  u tsp?  (0)  and 

(0}  coalesce.  The  fact  that  the  integration  around 
the  branch  cut  is  negligible  has  not  been  addressed  here 
due  to  length  constraints.  However,  the  accuracy  of  the 
asymptotic  results  with  numerical  results  docs  much  to 
support  this  claim. 

[5]  L.  B.  Fclscn,  "Radiation  from  a  uniax tally  anisotropic 
plasma  half-space,”  IEEE  Trans ,  Antennas  and 
Propagtion,  vol  AP- 1 1 ,  pp,  469-484,  1 963, 

[6]  L.  Felsen  and  N>  Marcuvitz,  Radiation  and  Scattering 
of  Waves.  Englewood  Cliffs,  NJ:  Prentice-Hall,  Inc., 
1973. 


